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Abstract. Gd3Ga50i2, (GGG), has an extraordinary magnetic phase diagram, where no 
long range order is found down to 25 mK despite Qcw ~2 K. However, long range order is 
induced by an applied field of around 1 T. Motivated by recent theoretical developments and 
the experimental results for a closely related hyperkagome system, we have performed neutron 
diffraction measurements on a single crystal sample of GGG in an applied magnetic field. The 
measurements reveal that the H — T phase diagram of GGG is much more complicated than 
previously assumed. The application of an external field at low T results in an intensity change 
for most of the magnetic peaks which can be divided into three distinct sets: ferromagnetic, 
commensurate antiferromagnetic, and incommensurate antiferromagnetic. The ferromagnetic 
peaks (e.g. (112), (440) and (220)) have intensities that increase with the field and saturate 
at high field. The antiferromagnetic reflections have intensities that grow in low fields, reach a 
maximum at an intermediate field (apart from the (002) peak which shows two local maxima) 
and then decrease and disappear above 2 T. These AFM peaks appear, disappear and reach 
maxima in different fields. We conclude that the competition between magnetic interactions and 
alternative ground states prevents GGG from ordering in zero field. It is, however, on the verge 
of ordering and an applied magnetic field can be used to crystallise ordered components. The 
range of ferromagnetic (FM) and antiferromagnetic (AFM) propagation vectors found refiects 
the complex frustration in GGG. 



1. Introduction 

Gadolinium Gallium Garnet, Gd3Ga50i2, (GGG), occupies a unique position among the 
geometrically frustrated magnetic systems. It has been a subject of intense research for the 
past three decades, as in terms of the geometry it lies somewhere between the well-studied case 
of stacked triangular lattices and lattices with a much more severe degree of frustration, such 
as the kagome lattice antiferromagnets and the pyrochlores. In GGG the magnetic Gd ions are 
located on two interpenetrating corner-sharing triangular sublattices, as shown in Fig. 1. It 
has been known for more than twenty years that GGG has an extraordinary low-temperature 
phase diagram [1, 2]. No long range magnetic order (LRO) has been found in GGG down to 
25 mK, LRO can be induced only by an applied magnetic field above ~ 1 T. The magnetic phase 
diagram depends on the external field orientation, as the critical temperatures and critical fields 



Figure 1. Positions of the magnetic Gd ions 
in a garnet structure. There are 24 magnetic 
ions per unit ceh, they are divided into two 
interpenetrating sublattices. 

for H\\ [100] and H\\ [111] differ by 30-40% [1]. Another pecuhar feature of the H — T phase 
diagram of GGG it that the phase boundary between the LRO and spin-hquid phases has a 
distinct minimum at T ~ 0.18 K, analogous to the minimum in the melting curve of ^He [3]. 
In this paper we report neutron diffraction measurements on a single crystal sample of GGG 
performed in an applied magnetic field at low temperatures. These measurements reveal that 
the H — T phase diagram of GGG has a much more complicated nature than previously assumed. 

The interest in Gd3Ga50i2 has been renewed recently by several developments. Firstly, 
a new mean-field approach to the theoretical studies of GGG has been developed for zero 
magnetic field [4, 5]. The theory suggests that at low temperature GGG is on the verge of 
achieving true LRO [6]. Secondly, the latest experimental and theoretical studies of a closely 
related "hyperkagome" system, Na4lr308, have considered the role of classical and quantum 
fluctuations in lifting the macroscopic degeneracy of the apparent spin-liquid state [7]. The 
similarity between the lattice structures of Na4lr308 and Gd3Ga50i2 is obvious and further 
comparison of their magnetic properties would be valuable. Finally, the reasons for GGG being 
a very efficient material for low temperature magnetic refrigeration have been reconsidered in 
view of the presence of a macroscopic number of soft modes in frustrated magnets below the 
saturation field [8]. 

2. Experimental details 

Polycrystalline GGG was synthesised from stoichiometric quantities of Gd203 and Ga203 by a 
solid diffusion reaction at T = 1400° C for 12 hours with intermediate regrinding. The procedure 
of regrinding was repeated until an X-ray diffraction study showed no impurity phases in the final 
product. In order to minimise neutron absorption we have prepared a sample with 99.98% of 
the non-absorbing isotope, ^^''Gd; this isotope was supplied by Oak Ridge National Laboratory. 
This polycrystalline sample was used previously for neutron diffraction measurements in zero 
field [9] and in an applied magnetic field [10]. A single crystal sample of GGG was then grown 
by the floating zone method using a two mirror infra-red image furnace. The low-temperature 
specific heat measurements performed on the single crystal sample [11] have shown results that 
are very similar to the results reported by Schiffer et al. [12], which implies an absence of a 
significant amount of impurity. 

The single crystal neutron scattering measurements were carried out using the DIO instrument 
at Institut Laue-Langevin, France. An 80 x 80 mm^ two-dimensional microstrip detector was 
used in the diffraction configuration. Typical intensities of the main nuclear Bragg peaks for 
the incident wavelength neutrons of A = 2.36 A were several hundred counts per second. We 
have measured the magnetic diffraction patterns at temperatures between 50 mK and 0.7 K in 





Figure 2. Field dependence of the integrated 
intensity of the three peaks with a strong 
FM component, (112), (220) and (440) as 
measured in a GGG single crystal for ff||[110]. 
The arrows indicate the fields of 0.82 T and 
1.65 T where the peaks exhibit a change in a 
rate with which the intensity rises with field. 
The intensity of the field-independent nuclear 
peak (004) is also shown for reference. 



Figure 3. Field dependence of the intensity 
of the two AFM peaks measured at different 
temperatures in a GGG single crystal for H\\ 
[110]. Top panel - (002) peak, bottom panel - 
an incommensurate peak nominally observed 
at ICi = (1.18,1.18,1.32) slightly off-centre 
of the 2D position sensitive detector. The 
peak is indexed as (1.35, 1.05, 1.37) when off- 
centring is taken into account. 



fields of up to 2.5 T. The horizontal scattering plane contained the {hhl) reflections so that an 
external magnetic field provided by a vertical cryomagnet was applied along the [110] direction. 
The results of the measurements in the (hkO) scattering plane for a field applied along the [100] 
direction performed on the same sample at Berlin Neutron Scattering Center, Germany have 
been reported previously [13]. 

Apart from describing our new single crystal results, we have also reanalysed previously 
reported neutron powder diffraction data [10] and indexed the field-induced incommensurate 
magnetic Bragg peaks. 

3. Results and Discussion 

Prior to the application of a magnetic field, we have attempted to find the relatively sharp 
low-intensity magnetic Bragg peaks, which have been detected in a powder diffraction pattern 
at temperatures below 140 mK at incommensurate positions, such as Qi = 0.64 and 
Q2 = 0.85 and several others. Recent theoretical results [4, 5, 6] advocated a scenario that 
in the zero field regime GGG upon cooling is very close to developing conventional long range 
magnetic order, which corresponds to the appearance of the sharp but not resolution limited 
diffraction peaks. These peaks give a correlation length of approximately 100 A, suggesting 
extended short range order (over 8 cubic unit cells of 12.3 A each). Our attempts to find peaks 
of the (hhl) type with a scattering vector length equal to Qi or Q2 in single crystal experiment 
were unsuccessful. 

The application of an external magnetic field at low temperature results in an intensity 
change for most of the peaks, which could be divided into the three main sets, ferromagnetic, 
commensurate antiferromagnetic and incommensurate antiferromagnetic. The ferromagnetic 
peaks have intensities that increase with the field and saturate at some higher field, while the 
antiferromagnetic peaks have intensity that grows in low fields, reaching a maximum at some 
intermediate field then decreases and (almost) disappears at fields above 2 T. 

Fig. 2 shows the field dependence of the integrated intensity, /, of the three FM peaks (112), 



Table 1. A description of the various k- vectors present in GGG powder profiles [10] at 
T = 0.08 K in different fields determined by reverse-Monte Carlo indexing [14, 15]. 



Field, T ki = (000) k2 = (001) = (0, 0, 0.724) 

1.0 • • • 

1.6 • • 

3.3 • - 



(220) and (440) as well as of one purely nuclear peak (004) measured at the base temperature 
of 50 mK. One immediate observation from this figure is that there are two fields at which the 
gradient of the I{H) curves exhibit a change: a pronounced flattening at 1.65 T presumably 
associated with magnetisation saturation and another change at half this field where the gradient 
increases by approximately 50%. The changes are most noticeable in the field dependence of the 
(112) peak, as the experimental point density there is the highest, but the same tendency 
is clearly visible for the (220) and (440) peaks as well. The data in Fig. 2 also indicate 
that there is no significant temperature dependence of the intensity of the (112) peak for 
50 mK < T < 250 mK. Nuclear contributions to the (112) and (220) reflections are allowed 
by symmetry, but this component is small compared to the FM field-induced intensity. 

Fig. 3 shows the evolution of the intensity of two AFM peaks (commensurate and 
incommensurate) with an applied magnetic field. These peaks are purely magnetic as the 
translational symmetry of the crystal structure prevents any nuclear contribution. One of 
the peaks, (002), clearly shows two local maxima at the base temperature of 50 mK. The 
positions of these maxima (0.8 T and 1.65 T) match well the fields where the I{H) curves 
for the FM peaks show a change in gradient. Another peak, ICi, is incommensurate, indexed 
as (1.35,1.05,1.37), but is nominally observed at a point (1.18,1.18,1.32) slightly off-centre 
on the 2D position sensitive detector. Unlike the FM peaks, the ICi is sensitive to small 
temperature variations (see Fig. 3 bottom panel). Apart from the AFM peaks shown in Fig. 3, 




Figure 4. Observed (symbols), calculated 
(full lines) and difference neutron powder 
diffraction profiles for GGG [10] in an applied 
field of (a) 3.3 T, (b) 1.6 T and (c) 1.0 T at 
r = 80 mK. 



Figure 5. Field dependence of the AFM 
peaks intensity, as measured by ramping an 
external field -ff||[110] up/down (solid/open 
symbols) with a rate of ~ 0.1 T/min. ICi^2,3 
denote the incommensurate peaks indexed in 
the main text. 




Figure 6. T dependence of the peaks 
intensity for i^||[110] at 1.24 T and 1.67 T (top 
and bottom panels). The arrow on the upper 
panel indicate samples cooling or warming. 
For the lower panel the data were collected in 
a field warming regime. The notations ICi^4 
correspond to the incommensurate peaks, 
their indices are given in the main text. 



Figure 7. A comparison of the magnetic 
phase diagram of GGG determined by bulk 
properties measurements (lines) [1, 2] and 
neutron difi'raction results (symbols). Solid 
and dashed lines correspond to the suggested 
regions of LRO and SRO. Apart from the new 
ILL results, neutron data from the previous 
HMI experiment [13] are also shown. 



we have observed similar behaviour for the half-integer peak (22^) ^.nd the incommensurate 
peaks IC2 = (-1.53,-1.53,0.51), IC3 = (0.59,0.59,1.95) and IC4 = (1.46,1.46,-0.51) indexed 
as (-1.48, -1.56,0.51), (0.52,0.68,2.04) and (1.41, 1.60, -0.51) respectively when off-centring is 
taken into account. The relatively small number of the magnetic reflections observed (largely 
limited to an [hhl) scattering plane) does not allow for a full structure refinement. However, 
even if a larger number of peaks was available, the magnetic structure determination would 
not be straightforward, as the AFM peaks appear and disappear in slightly different magnetic 
fields. This is well demonstrated in Figure 4 where the incommensurate and then commensurate 
antiferromagnetic peaks disappear upon increasing the field. The refined propagation vectors 
and their presence in the different field regimes are shown in Table 1. The peaks are not 
resolution limited, indicating that the field-induced magnetic order in GGG is short range in 
nature. 

The data in Figs. 2 and 3 were collected by fixing an external field and performing tu-scans 
for each refiection. This method ensures an accurate measurement of the intensity and the 
position of the Bragg peaks. Alternatively, a much more rapid data collection is achieved by 
constantly ramping an external magnetic field and simply counting the number of scattered 
neutrons per second. Although this method is less accurate in terms of measuring absolute 
peak intensity (particularly if the peaks are changing their positions or widths with field), it 
allows the study of any sample history dependence and observation of possible metastable states. 
The results of such an approach are summarised in Fig. 4. As can be seen from this Figure, 
the incommensurate peaks are very sensitive to the sample's field and temperature history and 
demonstrate significant hysteresis. Less pronounced hysteresis is also observed for a (002) peak 
between 1.1 T and 1.4 T. Remarkably, the (|^2) seems to exist only in the vicinity of (libO.2) T, 
while the incommensurate peaks have non-zero intensity over a much wider region between 0.8 T 
and 1.8 T. 

Having seen which peaks appear at base temperature in an applied magnetic field, it is 
equally important to follow the temperature dependence of these peaks. Fig. 5 shows the 
temperature dependence of the intensity of the selected FM and AFM peaks in a field of 1.24 T 



and 1.67 T. A large difference between the results for field-cooling and field- warming regimes 
for the incommensurate peaks again points to a metastable nature of the field-induced state. 
The instability and history-dependence were seen in the Monte Carlo simulations of GGG [16], 
although doubts were raised whether the energy of dipole-dipole interactions was calculated 
properly [4, 5]. 

Fig. 6 shows the magnetic phase diagram of GGG, where our neutron scattering results are 
compared to the bulk properties measurements [1, 2]. 

4. Conclusions 

To summarise, GGG is an intriguing magnetically frustrated system where it is difficult to 
project out a single ground state in an applied magnetic field. 
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